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Resumen. Is it possible to define, for certain values n the product of vectors of the real vector space of 
n dimensions , such that this is, with respect to multiplication and the ordinary addition of vectors, a 
numerical system which contains the system of real numbers? It can be proven that this cannot be done. 
CNJ ■ In the space of four dimensions this construction is possible if we are apart from the commutativity 

of the multiplication. The resulting system is the one of QUATERNIONS. In this work I first do 
a reminder of the fundamental concepts of Hamilton's Hypercomplex and then a deep work with such 
\ concepts. 
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1. Define quaternion 

Definition 1.1. We shall call quaternions or simply Hamilton's hypercomplexes to an expression in 
the form: 

^3 '. Q = a + bi + cj + dk 

where: a,b,c,d G KL Besides k are imaginary units, pairwise solutions of the equation x 2 = — 1, 
satisfying: 

ij = k = -ji 

jk = i = —kj 

ki = j = —ik 

ijk = -1 

i*=j* = k 3 = -l. 
CO ■ 

O ' Definition 1.2. We shall say a quaternion is purely imaginary, if the first element of the expression 

r-| is equal to zero (a = 0, Jm(H)). 
^— > ■ 

S Theorem 1.3. We say Q = Q' , with Q' — a' + b'i + c'j + d'k (Q,Q' G M), ie two quaternions are 
equal, if and only if, are equivalent the components of their imaginary and real parts: a = a',b = b',c = 
>' c\d = d'. 
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Proof. It's a direct consequence of the equality in M n . 

1.1. Fundamental Definition. 

• The sum and substraction are defined component by component, ie: 

Q + Q' = (a + a') + (b + b')i + (c + c')j + (d + d')k (1.1) 

Q-Q> =( a -a') + (b- b')i + (c - c')j + (d- d')k. (1.2) 

• The product is defined in the way: 

Q" = ( y a + bi + cj + dk){a! + b'i + c'j + d'k) (1.3) 

resulting: 

Q" = aa' - bb' - cd - dd! +(ab' + a'b + cd! -c'd)i + (ac' + a'c- bd' + b'd)j + (ad' + a'd + be' -b'c)k, 



with: 

a" = aa 



" ' bb' - cd - dd' 
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b" = ab' + a'b + cd' — cd 
c" = ad + a'c - bd' + b'd 
d" = ad' + ad + be — b'c 

Q" = a" + b"i + c"j + d"k. (1.4) 

1.2. Fundamental Properties. 

• Commutativity of the sum: Qi + Q 2 = Q2 + Qi- 

• Associativity of the sum: (Q 1 + Q 2 ) + Q3 = Qi + {Q2 + Q?,)- 

• Associativity of the product: (Qi • Q2) • Qz — Qi • (Q2 • Qz)- 

• Distributivity: (Q 1 + Q 2 ) ■ Q 3 = (Qi + Q3XQ2 + Qs)- 

Remark 1.4. Among the upper properties the following is missing: 

Q1Q2 ^ Q2Qi(ingeneral). (1.5) 
A very important result is Hamilton's hypercomplex are not a commutative field. 
Definition 1.5. We shall call conjugate of a quaternion Q, and denote Q, the number: 

Q = a — hi — cj — dk. (1.6) 

Let us write now the sum and the difference of a quaternion with its conjugate: 

• Sum: 

Q + Q = ( a + bi + cj + dk) + (a - bi - cj -dk) = 2a = 2® JRe(H). 

• Difference: 

Q -Q = (a + bi + cj + dk) - (a - bi - cj - dk) = 2(bi + cj + dk) = 2® Im(H). 

• Product by its conjugate: 

Q -Q = (a + bi + cj + dk)(a - bi - cj - dk) = a 2 + b 2 + c 2 + d 2 . (1.7) 

1.2.1. Other properties. 

• Selfpowered: 

• Additivity: 

Qi + Q2 = Q~i + Q~2 

• Multiplicativity: 

Qi • Q2 = Q2 ■ Qi 

• Divisibility: 

{ Q2 Q 2 

Is it possible to establish the inverse for the sum and the multiplication? 

Definition 1.6. From Q x + Q2 = we have: a + a' = 0, b + b' = 0, c + d = 0, d + d' = 0. This implies 

a' = —a, b' = —b, d = — c, d' = —d . Therefor Q 2 = —a — bi — cj — dk . Q 2 = —Qi- 

Analogously, 

Definition 1.7. If Q\ ■ Q 2 — 1 with Q x 7^ then Q 2 will be the multiplicative inverse for Q x . 



a-bi-cj - dk 
^ 2 a* + V + c 2 + d? 1 



with Qi = a + bi + cj + dk. 
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Now, we've already defined the conjugate and multiplicative inverse. Let's define division. 
Definition 1.8. Division is define like this: 

7T = Qi- Q2 1 (1-9) 

Remark 1.9. Notice when dividing we multiply the numerator by the multiplicative of the denomina- 
tor(recall (1.5)) 

2. Absolute Value 

Definition 2.1. We shall call absolute value or modulus of the quaternion, the nonnegative real 
number \Q\, 

\Q\ = \ a + bi + cj + dk\ = Va 2 + b 2 + c 2 + d 2 . 

Evidently, if we want to find the modulus of any quaternion: Q + Q' , Q", Q — Q' ,etc; this will be 
the square root of the sum of the squares of the real elements of each imaginary unit. Let us notice 
that: 

\Q\ 2 = Q-Q. 

As seen in formula (1.7). 

Thus, \Q\ 2 = a 2 + b 2 + c 2 + d 2 . 

Now with Q = a + bi + cj + dk and Q = a — bi — cj — dk (finding the conjugate of the earlier): 

• \Q\ = \Q\ 

• |Qi-Q 2 | = IQ1IIQ2I 

• IQiQ 2 | 2 = |Qi| 2 -|Q 2 | 2 . 

2.1. Norm. 

Definition 2.2. The norm will be defined like this: 

\\Q\\ 2 = Q-Q = \Q\ 2 (2.1) 

Now with Q = a + bi + cj + dk and Q = a — bi — cj — dk (finding the conjugate of the earlier): 

• IIQII = WQW 

• ||Qi-Q 2 || = IIQi||||Q 2 || 2 

• IIQiQ 2 || 2 = IIQi|| 2 -||Q 2 ||- 

It is simple consequence of the modulus. 
Remark 2.3. We defined division previously. Now we state the following equality: 

Qi _ Qi • Q2 
Q~2 ~ ||Q 2 || 

Definition 2.4. We shall call unit quaternion the Hamilton's hypercomplex which satisfies: 

IMI = 1 (2-2) 

i.e. a 2 + b 2 + c 2 + d 2 = 1. 
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3. Several ways of defining a quaternion 

3.1. Vector Form. 

Y : H — > M 4 
a + 6z + cj + d/c — >■ (a, 6, c, d); a, b, c, d G M 

Let us prove K is a bijective application. 
Vxi, :r 2 G EI with xi 7^ x 2 =>■ ^(^1) 7^ ^(^2)- 

- Therefor is injective. 

G M^3x G H : Y(x) = y. 

- Therefor is surjective. 

We've proved Y is a bijective function. 

Provided that(H, +) is an Abelian group, with the sum as defined earlier; that the product is 
distributive with respect to the sum, Q\(Q2 + Q3) = Qi • Q2 + Qi • Q3 and the product is associative 
Qi(Q2Qs) = (QiQ2)Q3', then, quaternions with the operations of sum and product define a ring. 

(H, +, *) is a ring. 

3.1.1. Definitions of this notation. 

(1) The quaternion(0, 0, 0, 0) is the neutral quaternion for the sum. It's obvious from the definition 
of sum: (a, b, c, d) + (0, 0, 0, 0) = (a, b, c, d). 

(2) Analogously, the neuter for the product is (1, 0, 0, 0) because (a, b, c, d)(l, 0, 0, 0) = (a, b, c, d). 

(3) Noting we are using our defined product: a(a, b, c, d) = (aa, ab, ac, ad). 

(4) (0, 1, 0, 0)(0, 1, 0, 0) = (—1, 0, 0, 0) as we can see this quaternion is identified by i 2 . 

From now on, we'll consider: Y(l) = (1,0,0,0) Y(i) = (0,1,0,0) Y(j) = (0,0,1,0) Y(k) = 
(0,0,0,1). 

We wonder, is it an isomorphism? We only have to prove the following holds: 

Y[(a + bi + cj + dk) + (s + ri + tj + hk)] = Y[a + bi + cj + dk] + Y[s + ri + tj + hk]. 

3.2. Other Vector Form. 

r : H — >■ 1R 4 
a + bi + cj + dk — >■ (a, v); a 6 8,7 vector in IR 3 



The conjugate of this vector is: (a, v ) — (a, — v ). I remark the multiplication with this notation 
is: Qi • Q2 — (a, ~v){o!, v±) = (aa' — v ■ v± , av\ + a' • v + v x Vi) 

3.3. Matrix form. The matrix form of defining a quaternion is: 



a + W + c ^' + dfc — ^(-c + di a-w) 



In order to illustrate this notation, it's convenient to develop the following, by taking 



«(i)=( i ?)' n(i) ^(o °)-^=(-i 



These are called Pauli's matrixes. 
The quaternion should be then: 



r^, , • ■ ,, \ (a + bi c + di\ 
n(a + ln + cj + dk)=(_ c + d . a _ u j 

fi(a + 6i + cj + dA;) = a- ^ ( l ) J^+^'^q -° 2 )+ c '(-°i o) +rf '(? 

Q(a + 6i + cj + dfc) = afi(l) + 6ft(i) + cQ(j) + dQ(k). 
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We have only left to wonder, is this function an isomorphism? It is, indeed, and the proof is equivalent 

to the earlier function (Y).li we find the determinant of the matrix f a + ^ c + ^ \ \ anc [ calculate 

v ' y —c + ai a — m J 

the modulus of the quaternion a + bi + cj + dk we can see identical results, therefor it is easily seen 

this function is an isomorphism. 

3.4. Trigonometric form. 

II : EI — > T ,T : Trigonometric form. 
a + bi + cj + dk — > pCisO + (p Cis(3)j. 

To illustrate this definition we write: 

a + bi = pCis9,c + di = p Cis/3,—c + di = —p Cis(—(3),a — bi = pCis{—6) with a + bi,c + di, —c + 
di, a — bi e C. 

Finding the modulus of the earlier complex number we obtain: 

p — \a + bi\ = \a — bi\ = \J a? + b 2 , p = \c + di\ = \ — c + di\ = \J c 2 + d 2 . 

3.5. New Trigonometric form. 

T : H — >■ J ,J :New trigonometric form. 
a + bi + cj + dk — >■ p(Cos(a) + Sin(a)j). 
a + bi + cj + dk — > pCjs(a). 

Remark 3.1. The following is the short way to express the earlier. We know the field of the quaternions 
is not commutative, so all of the transformations are equivalent. 

pCjs(a) = p{Cos{a) + Sin(a))j (3.1) 

The idea of the development is to shorten the trigonometric form. 

From complex analysis we know that: 

Ln(z) = Ln\z\ + iArg(z) 

Tarr\z) = -Ln( — ) 
2 i — z 

Definition 3.2. Applying Pythagorean Theorem in R 4 , ie, working with the complex axes(Fig.l), we 
define: 

p= ^{a + bi) 2 + {c + di) 2 (3.2) 

c ~\~ d% 

a + bi = pCos(a); c + di = pSen(a); Tan(a) = — — — (3.3) 
From (3.3) we have that: 

a = Tan~\ C -±4) = Cos~\ a -±^) = Sl n'\ C -±^) 
y a + bi p p 

Now, for obtaining the true value of (a) any of the three following trigonometric function value's must 
be calculated: Tan -1 , Cos -1 , Sin -1 . They must be equal to (a). Now we wonder, why not to consider 
the value of (p) as the modulus of the quaternion? The answer is the following: The trigonometric 
equality(3.10) does not always hold, only when: b = d = 0. And the hypercomplex were not defined 
under these conditions. The conclusion I reach about it, is that if these equalities hold then we can 
take the modulus of the quaternion as (p). From (3.7) we have that: 

m -i,c + di. i T ,a + d+(b~c)i. 

Tan — = -Ln{ — — 

a + bi J 2 a- d+ (b + c)i J 



6 



MIJAIL ANDRES SARALAIN FIGUEREDO 



ie: 

1 a + d + (b — c)i. 
a — T^Lni ; — — ) 

2 K a-d+ (b + c)i J 

As we already can see we have p which gives us (3.8), (a) which yields us (3.10) or (3.12). Then we 
can write a hypercomplex as: 

a + bi + cj + dk = pCjs(a) 

where p, a G C. 



In the earlier graph we can see that a quaternion can be represented in a system of two complex 
coordinates. It is impossible to represent the set C in a coordinate axis. I only do it to see where I 
obtain (3.8) and (3.9). At the same time, each complex axis I show has values in R 2 ,ie, we have the 
earlier graph in R 4 . 

The remark I will write further is just a question. When defining a field in trigonometric form, do 
we have to reference the before field? 

As we can see I've arrived to a new way of defining quaternion. Now we have: 

pCjs(a) G H (3.4) 

with p, a G C, then we can write p = p\{Cos((3) + iSin((3)) with p\ G R, a G C. So 

a + bi + cj + dk = pCjs(a) = pi(Cos(/3) + iSin{f3)){Cos{a) + Sin(a)j) 

From complex analysis we know that: 

Sen(a) = Sen(x + iy) = Sen(x)Cosh(y) + iCos(x)Senh(y) 

Cos(a) = Cos(x + iy) = Cos(x)Cosh(y) — iSen(x)Senh(y) 

Substituting 

a + bi + cj + dk = pi[Cos(/3) + iSin(/3)][Cos(x)Cosh(y) — iSen(x)Senh(y) + [Cos{x)Cosh{y) — 
iSen(x)Senh(y)]j]. 
yields: 

a = piCos(x)Cosh(y) 
b = —piSin(x)Sinh(y) 
c = piSin(x)Cosh(y) 
d = piCos(x)Sinh(y) 

William Rowan Hamilton in " On a New Species of Imaginary Quantities Connected with a Theory of 
Quaternions", wrote: 

a = pCos{9) 
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b = pSin(9)Cos($) 
c = P Sin(6)Sin(6)Cos{i)) 
d = pSin(9)Sin(e)Sin(^) 

3.6. Trigonometric Matrix Form. 

* : H — >■ M 

a + bi + d + dk^ ( P Cis9 P° Cis P 
a + bi + cj + dk^ ^ _ potts /3 pC7is(-e; 

I develop this kind of definition mostly for a practical work, such as multiplying or dividing, due to 
these operations are easily performed in trigonometric form. 
Let's see that: 

pCisO Po Cis[3 \ _ ( pCisd \ / p Cis(3 
-poCisp pCis(-e) J ~ \ pCis(-e) J + \ -poCisp 
pCisQ p Cisp \ _ f pCosO \ / ipSend 
-poCis[3 pCis(-6) J ~ \ pCos(-e) J + \ ipSen(-9) 
poCosp \ ( 
-p Cosf3 y V -poSenfi 

3.7. Logarithmic Form. Let us define the following map: 

r : H — > L 
o + 6i + cj + <i/c — > In ab l c>d k 

This way of definition is to convert the elements a,b, c, d G K in natural logarithms. 

a = In a, b — In 6, c = In c, = In d 

(the upper is an abuse of notation) Obviously it's possible to think each part of the quaternion as real 
logarithms, because they are real. This notation is only possible when the real elements are not zero. 

3.8. Exponential Form. 

A : H — >E 
a + bi + cj + dk — > pe %e + poe tf3 j. 

Let us define the exponential form . First, we know from complex that:e z = e x (Cosy + iSeny) ,this 
means that 

e yi _ Q 0S y _|_ iSeny = Cisy. (3.5) 

Then: CisB = e ie , Cis{3 = . 

As we know a a + bi + cj + dk can be written as pCisO + p Cis/3j. 

a + bi + cj + dk = pCisO + p Cis/3j, by (3.5), a + bi + cj + dk = pCisO + p Cis/3j = pe td + poe tl3 j. 
Now calculate (a + bi + cj + dk) n . 

(a + bi + cj + dk) n = (pe td + po^i^developing Newton's binomial we have: (pe td + poe^j)" 1 = 
= {pe ie ) n + Cl{pe ie ) n ' x p G e ifz 'j + Cl{pc ie ) n - 2 {p G e lli jf + ...CI' 1 pc ib \p G e if) 'j) n ~ l + (p e i/3 j) n - 

(pe W + Poe^T = E (l ) (P eid ) n ~ h (P^j) h (3-6) 
h=o ^ ' 

(pe* + poeVj)" = (pe W ) n £ ( £ ) (^e^j) h (3.7) 
Now we must verify by induction if (3.7) holds. 

Obviously it holds in the case n — 1. Suppose it holds for n = k and prove it is true in the case 
n — k + 1. Now, with the new definition in trigonometric form we'll prove the following: 
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Theorem 3.3. We know a+bi+cj+dk = pCjs(a), then (a+bi+cj+dk) n = (pCjs(a)) n = p n Cjs(na). 

4. Functions of a hypercomplex variable 

If a variable w is related with z such that to each value of z in EI corresponds a value or set of defined 
values of w, then w is a function of the hypercomplex variable z, w — f(z). If z = a + bi + cj + dk and 
w = u + v i + sj + tk with the values of a, b, c, d, u,v, s,t G R. u + v i + sj + tk = /(a + bi + cj + dk), 
and each of the real variables u, v, s, t e R are determined by the real quartet a, b, c, d. That is to say, 
u = u(a, b, c, d),v = v(a, b,c,d), s = s(a, b, c, d),t = t(a, b, c, d). 

Example 1: w = z 2 + 5. 

u + vi + sj + tk = (a + bi + cj + dk) 2 + 5, 

u + vi + sj + tk = a 2 - b 2 - c 2 - d 2 + 2abi + 2acj + 2adk + 5. Then: 

u(a, b, c, d) = a 2 — b 2 — c 2 — d 2 + 5, 

v(a, b, c, d) = 2a, s(a, b, c, d) = 2ac, t(a, b, c, d) = 2ad. 

4.1. Limit Definition. Let f(z) a function defined in all the points in some neighborhood of z . We 
say that Wo is the limit f(z), when z tends to zo, 

lim f(z) = w . 

Z-+ZQ 

That is, for all positive epsilon exists a positive number lambda such that: 

\f(z) - w \ < e 

when: 

\z - Zq\ < \{z ^ Z ). 

Suppose that, 

lim f(z) = u + v i + s j + t k 

z^zo 

where f(z) = u + vi + sj + tk, z = a + bi + cj + dk, z = a + b i + c Q j + d k. Then by the inequality 
it becomes in: 

\u + vi + sj + tk — (u + v i + s j + t k) \ < e 

when: 

| a + bi + cj + dk — (a + b i + c j + d k) \ < A. 
Making algebraic transformations, (4.5), (4.6) yield: 

\(u - u ) + (v- v )i + (s - s )j + (t- t )k\ < e 

when: 

| (a - a ) + (b - b )i + (c - c )j + (d - rf )/c| < A. 
But from (4.7) and (4.8) we have, applying (2.1): 

y/(u - u ) 2 + (v- v ) 2 + (s- s ) 2 + (t- t ) 2 < e 

when: 

V(a - a ) 2 + (6 - 6 ) 2 + (c - c )2 + (d - d ) 2 < A. 
Uniqueness of the limit of a hypercomplex function: Suppose there exist two limit points wo,wi 
(wo 7^ Wi). By definition of limit: 

\f(z) - w Q \ < |, t«/ien|2; - z \ < A 

|/(2) — wi\ < -, when\z — z x \ < A 
Let's work with the real number^! — w \. The following step is obvious: 
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\wi ~ W \ = \W! -W + f(x) - f(x)\ = | - f(x) + Wi + f(x) -W \ = \- (f(x) - U>i) + f(x) - Wq\ < 

e e 
2 + 2 



< | - (f(x) - Wl )\ + \f(x) - w \ = \ f(x) - Wl \ + \f(x) - w \ < J + I < e 



\wi — wq\ < e, 

but as u>o,wi are constant \wi — wo\ cannot be made as small as one wants. Then, wo = w\. 
Define a function f(z) = w = u + vi + sj + tk such that: 

s ,. u(a + Va) — u(a) ,. via + Va) — v (a) . ,. s(a + Va) — sia) . ,. t(a + Va) — i(a) , 
/'(z) = hm — — + hm — —i+ hm — — 7 + hm — —k. 

Va^o Va Va^o Va Va^o Va Va^o Va 

r// n ,. u{b + Vb)-u(b) u(6 + V6) - . s (6 + V6)-s(6) . ,. *(& + V6) - t(6) , 
f(z) = hm — — — + hm — — f — s+ hm — —f — 7 + hm —f —k. 

vb^o V6z V6^o Vbi V6^o Vbi V6^o Vbi 

= Um U(c + Vc)- U (c) + Um «,(c + Ve)-»(c), + Um .(c + Ve)-,( C ) ^ t (c+ Vc) - t (c) fc 

Vc-*0 Vcj Vc^O Vcj Vc->0 Vcj Vc^O Vcj 

= Um 1 ,(d + w)- MM+ Um ,(d + w) -, (<0 Um s( d + w) -.(d) Um ^ + w)-*W, 

vrf^o VaA; vd^o Va/c vd^o VaA; vd^o VdA; 

Then: 



du dv . ds . dt 
da da % da^ da 


(4.1) 


,,, s du . dv ds . dt . 
fiz) = -db l+ db + db k -db J - 


(4.2) 


f'( ^ U ' 9v ^ ^ ds ^ dt . 
dc^ dc dc dc 


(4.3) 


fl( s du dv . ds . dt 

I[z) = ~dd k+ dd J ~dd t+ dd- 
conclude that a hypercomplex function is 


(4.4) 

analytic or integer if: 


du dv ds dt 
da db dc dd 


(4.5) 


dv du dt ds 
da db dc dd 


(4.6) 



ds dt du dv 
da db dc dd 

dt ds dv du 
da db dc dd 

We call the latter the generalized theorem of Cauchy-Riemann. 

References: 

[1] Kurochov, Algebra Lineal. 

[2] Teresita Noriega, Algebra Lineal II. 

[3] Ebbinghaus et al, Number, Stringer, 1991. 



(4.7) 
(4.8) 



10 



MIJAIL ANDRES SARALAIN FIGUEREDO 



[4] Peter Grogono, Rotation with Quaternions, December 2001. 

FACULTAD DE MATEMdTICA, FlsiCA Y COMPUTACIoN, UNIVERSIDAD CENTRAL MARTA AbREU DE LAS VIL- 
LAS, APARTADO Postal 54830, Santa Clara , Villa Clara, Cuba 
E-mail address: mi jail@mf cuclv.edu. cu 



